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EINSTEIN OBSERVATORY RESULTS ON ACTIVE GALACTIC NUCLEI* 


R. F, Mushotzky and S. S. Holt 

NASA/Goddard Space Flight Center, Greenbelt, MD 20771, U.S.A, 


ABSTRACT 

The Einstein Observatory, by virtue of its increased sensitivity and improved 
angular resolution, has increased substantially the number of known X-ray emitting 
active galactic nuclei. This has made possible the detailed study of the relation 
of X-ray flux to both the continuum flux in the optical and radio bands, and to the 
line emission properties of these objects. In addition, the Einstein imaging 
instruments have detected morphology in AGN X-ray emission, in particular from 
jetlike structures in Cen-A, M87, and 3C273, 

The improved energy resolution and sensitivity of the spectrometers onboard the 
Observatory has provided detailed information on the geometry and ionization 
structure of the region responsible for the broad optical emission lines in a few 
AGN's. Such information, combined with detailed theoretical model i no and IUE and 
optical observations, have allowed us to construct a moderately detailed picture of 
the broad line region in these objects. 

INTRODUCTION 

The Einstein Observatory consisted of a focusing X-ray telescope with four focal 
plane {and one non-focal plane) instruments [1]. The two imaging focal plane 
instruments were the High Resolution Imager (HRI) with ~ 4" resolution and ~ 10 cm 2 
sensitive area in roughly the .2-2 keV band, and the Imaging Proportional Counter 
( IPC) with ~ 60" resolution and 80 cm 2 sensitive area in the .2-3.5 keV band. The 
focal plane spectrometer with the most extensive results on active galaxies is the 

^Review of literature finished April 1982 
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Solid State Spectrometer (S$S) which had a 6‘ diameter beam, 160 eV FWHM energy 
resolution, and ~ 100 cm2 of effective area in the .5-4.0 keV bend. 

Because of the low background per spatial resolution element, made possible by the 
focusing optics, the Einstein Observatory had a threshold ~ 100 times lower than 
that of the HEAO-1 satellite for source detection by the Imaging detectors and a 
threshold for detailed spectroscopy for sources at about the HEA0-1 limit. However, 
because most previous X-ray experiments were most sensitive in the 2-20 keV band, 
the overlap in energy band between Einstein and previous X-ray surveys is small. 

During its lifetime, the Einstein Observatory looked at over 5000 targets of which ~ 
1000 were active galaxies. In addition to the targets which were selected for 
observations, there were "everal hundred "serendipitous" observations of active 
galaxies that were in the ~ 1 square degree fields selected for other reasons. 

SURVEYS 

The main result of the Einstein observations of AGN, in our opinion, has been the 
confirmation of the ubiquity of X-ray emission from these objects; that is, X-ray 
emission is now established as a basic characteristic of AGM. This result, combined 
with the broad band HEAO-1 spectra, shows that much of the luminosity of the AGN 
often appears in the X-ray band. In fact, the several hundred active galaxies 
discovered because they were X-ray sources in Einstein imaging fields is a 
significant contribution to the total number of known AGN. 

The vastly increased data base on X-ray luminosities from Einstein observations has 
enabled a number of optical, radio, and X-ray correlations to become well 
established. While several of these correlations were hinted at in pre-Einstein 
data [2,31 the large dispersions in the correlations require the larger HEAO-2 data 
base to confirm them. 

A. X-ray and Optical Luminosities . In two large survey papers, Zamorani et al. 
(1981) [4] and Ku, Helfand, and Lucy (1980) [5] (KHL), it was shown that for 
optically selected AGN with m y < 18.5 and log f > -12.2 erq/cm 2 sec the optical and 
X-ray fluxes are strongly correlated. While there is a fair dispersion in the 
relation, the .5-3.5 keV X-ray luminosity L is typically - .2 of the optical 
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luminosity L 0 for an AGN. The values range from upper limits of L x ~ 3xlO“'\ Q to a 
maximum of L^ ~ 3L q . Since the 3.5-60 keV luminosity of an AGN is usually 3.5 times 
its .5-3.5 keV luminosity, the typical intrinsic X-ray and optical luminosities are 
similar. 

Recent work by Zamorani [61 and Reichert et al, [7] show that the X-ra.v to optical 
luminosity ratio decreases as a function of optical luminosity, with the higher 
luminosity objects havinq a lower X-ray to optical ratio such that L x a l 0 * 7 . This 
relation appears to hold over roughly five orders of magnitude, between low 
luminosity AGN (Se.yfert galaxies) at one end and QSO's at the other, and thus may 

argue for a continuity of properties between these objects. The L x /L 0 re1atl0n 1S 
an important constraint on theories of the enerqv generation mechanism in AGN's. 

Such theories must explain why QSO's have about an order of magnitude less X-rav 
emission per unit optical emission than Se.yfert galaxies despite their strong 
spectral similarities. 

B. X-Ray and Radio Properties . Many of the AGN in the initial surveys were radio 
loud objects, i.e., l n ~ 2 x 10" 3 L . Several authors have found that the X-ray to 
optical ratio is related to the radio flux such that radio briqht quasars have about 
triple the average L x to L 0 ratio, with Lx/L 0 als0 bein ^ a function o f the radio to 
optical luminosity Lr/l 0 . Using the objects in the KHL and Zamorani et al. surveys 
Setti and Woitjer [8] have shown that for optically bright AGN's (L 0 > 5 x 10^ 
ergs/sec) L x /l 0 increases roughly linearly with Lr/Lq. F. Owen and co-workers 
[9,10] have shown that for objects selected on the basis of their strong hi qh 

frequency radio flux, L r /l x shows a very small dispersion. Since the radio emission 
is due to synchrotron radiation from relativistic particles these correlations argue 
for a component of the X-ray radiation being due to emission from relativistic 
particles. These simple facts are complicated by the fact that the Setti -Woitjer 
relationship does not hold for low luminosity AGN's and that the Gwen relation does 
not hold for a sample of optical or X-ray selected quasars. In addition, the 
variance in the Setti -Woitjer relation is quite large. We believe that this argues 
for a non-simple relationship between X-ray emission and relativistic particle 
population, KHL [5] argue that synchrotron radiation seems to be the dominant 
emission mechanism from the IR to X-ray for nearly all h i qh luminosity AGN, with the 
radio bright objects having an additional X-ray emission component due to inverse 
Compton processes, Detailed X-ray spectra of BL Lac objects from the Einstein SSS 
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and HEAO-1 (a subset of radio briqht AGN ) tend to confirm the synchrotron origin of 
.2-4 keV X-rays from these objects. However, at present we do not have X-ray 
spectra of even a small sample of radio loud hi qh luminosity AGN's. 

G. X-Ray and Optical -UV Line Emission. It was noticed by Elvis et al. [21 based on 
the relatively small sample of AGN available in 1978, that some optical emission 
line properties are correlated with the X-ray luminosity. With the increased 
Einstein data base, several of these correlations are now well established. For 
example, Grindlay et al. [11], Kriss et al. [12], and Steiner [13] have shown that 
both the absolute He luminosity and the OIII/Hp ratio are well correlated with L x 
for radio quiet objects. C. C. Wu and co-workers [14,15] have noted that the CIV 
luminosity is strongly correlated with the X-ray luminosity. Since these parameters 
are also well correlated with optical luminosity it is not surprising that they 
correlate with X-ray luminosity as well [7]. 

Detailed photoionization models of the broad line region in AGN have been able to 
qualitatively reproduce the correlation of the permitted lines with the optical and 
X-ray continuum levels. Because many of the low luminosity AGN have extremely 
similar X-ray and UV conti nua (Wu [15], Mushotzky [16]), it is possible to treat 
them as a homogenous qroup. However, it is not clear at present how to produce 
either the correlations of 01 1 1 /He with L x j n pel I strong objects [13] or the lack 
of such correlation in Fell weak objects. Perhaps quantities other than the X-ray 
to optical ratio or the local spectral slopes are important, such as geometry or 
filling factor [17], or that selection effects are responsible for the apparent 
correlation. 

The positive correlation of the infrared to soft X-ray ratio with the X-ray 
luminosity suggests that the IR continuum is influenced by the X-ray flux. Lawrence 
and Elvis [17] argue that there is a short wavelength non-thermal IR component that 
is directly correlated with the X-ray flux, and a long wavelength IR component due 
to dust heated by a source indirectly related to the X-ray emission. 

Kriss et al. [12] have shown that there may exist a weak correlation of optical line 
width to X-ray flux; however, this depends mainly on including Seyfert II 1 s in the 
sample. Since absorption by dust and qas can have a strong effect on the width of 
He lines as well as on the Einstein fluxes (for it has been shown [16,17] that L x 
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and Nfl are anti -correlated) , we consider this relation suspect, and it is most 
likely due to the fact that the broad wings of He can be strongly reduced by 
reddening in the low luminosity objects. This observational selection effect 
probably also causes part of the OIII/Hu versus L x correlation since observed H 3 
would be reduced in even slightly reddened objects. For L K > 5 x 10 43 erg/s, the 
roughly constant width of Hp indicates that the mass motion of the line emitting gas 
is not directly related to the production of X-rays. 

JETS 

The X-ray emission from almost all active galaxies is not resolved by X-ray 
observations since the source of the radiation is probably smaller than a light year 
in size (3 x 10" arc sec at 20 Mpc), but deep exposures with the HRI and IPC have 
revealed X-ray emission from jets in Cen-A, 3C273, and M87 [18,19,201. While the 
jets in M87 and 3C273 were well known before the Einstein observations, that in Cen- 
A was discovered in the X-ray band [22], In all three cases the X-ray, optical, and 
radio jets have more or less the same structure [21], and the ratio of X-ray to 
optical fluxes is consistent with a synchrotron origin for the X-ray flux. However, 
this requires a steepening in the IR-optical region of the spectrum and thus 
Indicates that the relativistic electron population of the jet likewise steepens at 
high energies. This is similar to the radio-optical -X-ray spectra of BL Lac objects 
and is consistent with the conjecture that BL Lac objects are jets seen end on [23]. 

The mere presence of jet related X-ray emission is strong evidence for 
reacceleration of particles in the jet with relativistic factors ~ 107 . The 
detailed structure of the jet in Cen-A [22] (the only gne resolved by Einstein) 
independently argues for in situ acceleration (or reacceleration) since the size of 
knots ~ 10 arc sec (or ~ 250 pc in Cen-A) is so large that the X-ray emitting 
electrons could not live to cross the knot. The situation is even worse for 
relativistic electrons to cross from knot to knot or from the nucleus to the closest 
knot (a distance of ~ 750 pc). Evidence of hot gas in the central region of Cen-A 
[22] is roughly consistent with confinement of the jet by the thermal 
pressure of a hot (T ~ 107K) X-ray emitting plasma. Detailed X-ray spatial and 
spectral observations may give us information on the mechanisms of particle 
acceleration and transport and the mechanism of jet confinement. 
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Analysis of simultaneous IUE and SSS spectra of BL Lac obiects show that, in at 
least three cases, the X-ray and UV spectral indices are similar and that extension 
of the UV spectrum can roughly predict the X-ray flux C23]. These data have been 
interpreted as a synchrotron emitting relativistic Jet producing both the UV and X- 
ray radiation similar to the resolved jets of M87 and Cen-A. 

"NEW" TYPES OF ACTIVE GALAXIES 

Because of the ubiquity of X-ray emission from AGN, it is possible that "new 1, active 
qalaxies can be discovered by their X-ray emission even thouqh they may not show 
obvious optical or radio peculiarities. The optically "dull" galaxy NGC 4156 [24] 
may be such an object. Despite its proximity to the famous Seyfert NGC 4151, it was 
not discovered to have a active nucleus until Einstein's IPC and HRI observations 
showed a pointlike X-ray source of L x ~ I0 42 ergs/sec which is - .1 of the X-ray 
luminosity of NGC 4151. This object has very weak optical emission lines and a 
"normal" continuum. This discovery shows the weakness of any individual survey 
method for AGN and the power of observing in a new frequency range in an unbiased 
manner. The "serendipitous" Einstein sources give us the opportunity to discover 
objects or properties of objects which we have not been clever enough to 
anticipate. In fact, the recent history of astronomy has been largely one of 
serendipitous discoveries (X-ray and radio pulsars, quasars, X-ray emission from 
clusters of galaxies, etc.) 

The Einstein data on the very weak active nucleus in the "normal" qalaxy MR1 [25] 
corroborates the radio and optical data. The X-ray spectrum of MR7 [26], in 
particular, requires a component in addition to the thermal qas which is 
indistinguishable in form from that of more luminous AGN (see below). This 
establishes a connection between QSO's and the nuclei of some "normal" galaxies over 
a range of 107 in luminosity for the active galaxy phenomena. It is of interest 
that the lowest luminosity AGN X-ray source is only ~ 10 times as luminous as the X- 
ray binary SMC X-l* 

A serendipitous Einstein observation has revealed the first "narrow line" quasar. 
Stocke et al. [27] report that this X-ray source has a reddish continuum and 
relatively narrow 600 km/sec permitted lines. They argue that this object could be 
a prototype for a "new" class of "quasars" analogous to high luminosity Seyfert 2's. 
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It Is clear from this very limited list that X-ray observations, whatever their 
biases, can find active nuclei in over a wide range of types and luminosities. 

DETAILED SPECTRAL STUDIES 

Because spectrometers are, in general, less sensitive than imaqinq experiments, the 
Einstein Observatory spectrometers were restricted to looking at the relatively hi qh 
flux sources discovered by previous X-ray satellites. As reported by Holt r 281 
solid state spectrometer observations of active qalaxips confirm that the a = 0.65 
energy index of the Seyfert I galaxies in the 2-60 keV band established by Mushotzky 
et al. [29] from HEAO-1 data extends down to 0,6 keV. Thus a universal power law 
form for AGN appears to hold over at least 2 decades in enerqy from 0.6 to 60 keV, 

The most detailed results have been published on NGC 4151 [30]. The SSS data for 
this source do not show the form expected from a power law X-ray source extinguished 
by a uniform absorber: there is too much X-ray flux in the .5-2.0-keV band. The 
data can be well fit by a model in which the absorber is the broadline clouds, but 
their coveraqe is patchy so that X-rays can pass unextir.quished through holes in the 
cloud cover. Since the ratio of 7.1 keV Fe absorption to 6.4 keV Fe fluorescent 
reemission suggests spherical symmetry in the absorbing material, the model fit to 
the data determine a mean cover! no factor of the clouds of ~ .9. This large 
covering factor is somewhat surprising in view of the .1 coverinq factor inferred 
from Lyman continuum measurements for high redshift quasars. In addition, 
determination of the energy of the Si, S, and Fe {from HEAO-1 data) K-edqes 
demonstrate that the absorbing gas is cold. 

This model essentially established the global geometry of the clouds, titJfi r 
approximate sizes, and their ionization state. It has also proven to be a key to 
understanding that the anti -correlation of X-ray absorption and luminosity in 
Seyferts discovered by HEAO-1, and possibly the anti -correlation of CIV equivalent 
width versus luminosity (the "Baldwin" effect) in hT qh luminosity AGN's, are 
probably covering factor phenomena. The hiqher luminosity sources have, on the 
average, lower covering factors and, then-fore, less X-ray absorption and weaker UV 
1 i nes . 
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Detailed analysis of more SSS data Is proceeding and at least one other low 
luminosity galaxy has a large covering factor. The upper limits on the column 
density for high luminosity sources from the SSS is consistent with the coverlna 
factor being Inversely correlated with Intrinsic luminosity. The SSS data do not 
exhibit X-ray emission lines from AGN; this lack of evidence for oas at temperatures 
between 106.4 and 107.8 °k is consistent with recent models of the confining medium 
In the region emitting the broad optical emission lines. 

THE COSMIC X-RAY BACKGROUND AND AGN 

It was hoped that the vastly improved sensitivity of the Einstein Observatory would 
"solve" the "problem" of the Cosmic X-ray Background (CXB) either by resol vinq It 
into point sources or by establishing the contribution of different classes of 
objects to it. 

Well before the launch of Einstein, it had been suggested [311 that AGN were most 
likely to be the dominant point source contributor to the CXB. From the Initial 
Einstein survey of optically selected QSO's, it seemed, indeed, that an evolving 
population of AGN's (quasars) could easily make IDO percent of the ~ 1 keV CXB. In 
fact, as pointed out by Setti and Woltjer [323 and others, the CXB data seemed to 
require a turnover in the optical quasar counts at m ~ 20 in order that the CXB not 
be overproduced. 

A variety of complications have conspired to prevent the unambiguous resolution of 
the CXB, however. First, the spectral data from HEAO-l [33] have established that 
the CXB spectrum is much flatter (a ~ .4) than those of Seyferts between a few keV 
and a few tens of keV; deZotti et al. [34] have shown that the sample of Seyfert 
galaxies from which we have data cannot be representative of the sample which might 
produce the XRB, even when redshifted to z > 1. 

In addition, further work on quasar counts [35] and the radio-optical -X-ray 
correlations discussed in Section II showed that the problem has added complexity. 
While the last word is by no means in, the latest work [6,34,36,37] seems to 
indicate that the strongly evolving population of luminous (M < - 23.8) quasars make 
up less than 45 percent of the CXB with most likely values nearer 20 percent 
(however different authors disagree somewhat such that the calculated contribution 
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of all AGN to the CXB could be from CO to 10 percent). A point which is usually 
glossed over is that the Marshall et al CXB spectrum [33], consistent with a 45 keV 
bremsstrahlung spectrum, has the preponderance of the CXB energy density more than 
an order of magnitude higher In energy than the Einstein measurement. We miqht 
construct a scenario, therefore, whereby early AGN with spectra steeper than that 
typical of Seyferts might dominate the CXB near 1 keV and yet contribute virtually 
nothing to the total CXB energy density [38]. 

In order to truly resolve the problem one needs the spectra of AGN, particularly at 
energies > 10 keV, their luminosity function as a function of epoch, and better 
optical counts. Thus, despite the larqe Increase In Information on AGN provided by 
the Einstein Observatory, an unambiguous explanation for the CXB still remains a 
challenge 
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X-RAY SPECTRA OF YOUNG SUPERNOVA REMNANTS 


$.$. Holt 

Laboratory for HI qh Energy Astrophysics 
MASA/Goddard Space Fliqht Center 
Greenbelt, Maryland 20771 


A. HISTORICAL NOTES 

The early identification of the strong X-ray source in Taurus with 
the Crab nebula (Bowyer et al. 1964) was the first milestone in the 
association of X-ray emission with supernova remnants. Unfortunately, 
it proved to be "red herring" which clouded the interpretation of X-ray 
emission from supernova remnants for a decade. Because the Crab was one 
of the brightest X-ray sources in the sky at a few keV, the 
Interrogation energy of the early surveys, and because it was the first 
(and for several years the only) X-ray source conclusively identified, 
the potential association of a supernova origin with the large body of 
unidentified X-ray sources was not an unreasonable hypothesis. 

By 1970 proportional counter spectra from rocket-borne 
observations, with resolving power E/aE < 2, were available for Cas A 
and Tycho as well as the Crab (Gorenstein, Kellogg and Gursky 1970). 
Since the statistical significance of the data points were 
exposure-limited, they were fittable with either bremsstrahlung conti nua 
with temperatures of ~ 107 K or power-law spectra; the analogy with the 
Crab suggested to many interpreters of the data that the latter 
representation of the spectra were more appropriate. By that time the 
Crab spectrum s*a$ well -measured to be a power-law With a ~ 1 between 
1-100 keV, and relatively welh-understood (from the point of view of 
available enerqy, if not energy transfer) as arising from the 
non-thermal conversion of rotational enerqy from its pulsar. Similar 
pulsars which were unobservable because of their beaming characteristics 
were suggested for other young supernova remnants, since only youna 
pulsars could provide the required rotational energy loss (e.g. Pacini 
1971). Similarly, models featuring hidden pulsars In sources not 
associated with traditional supernova remnants (e.g. Davidson, Pacini 
and Sal peter 1971) were relatively popular. 

The discovery of "X-ray pulsars" (Giacconi et al. 1971), 
distinguishable from radio pulsars like the Crab by their inability to 
account for the source luminosity from rotational energy loss, provided 
the key to understanding the energy generation mechanism of the non-SMR 
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sources in terms of mass accretion onto neutron star components of 
binary systems, as had been suggested by Shklovsky (1967), But what 
about the young SNRs Cas A and Tycho? An early suggestion by Heiles 
(1964), that the SN blast wave sweepino up the interstellar medium into 
a shell of shocked material would provide an observable X-ray source, 
seemed more appropriate to older remnants in which enough mass had been 
swept up to provide an adequate emission measure, and enouqh time had 
elapsed to allow the radiating shell to equilibrate at a temperature low 
enough to cool efficiently (e.g, Woltjer 1972). Nevertheless, a blast 
wave interpretation for the young SNRs increased in popularity as the 
newer data indicated better agreement with a thermal spectrum. The key 
experimental result was the discovery of Fe line emission from Cas A 
(Serlemitsos et al. 1973) which, although originally interpreted by its 
discoverers as arising from charge exchange with ~ 10 MeV/nucleon cosmic 
rays, was quickly recoqnized as the signature of thermal processes in 
the remnant. As data from satellite experiments became available 
(Charles et al. 1975$ Pravdo et al. 1976; Davison, Culhane and Mitchell 
1976) this picture was refined to require two separate thermal 
components for consistency with these higher statistical precision 
data. Itoh (1977) suggested that the harder of these two components was 
associated with the blast wave, while the softer was associated with a 
thin, dense shell of ejecta accumulated from the "reverse shock" 
postulated by McKee (1974). The history of Cas A and Tycho measurements 
in the decade preceding the first spectral data from Einstein is given 
in Table 1. 


Table 1 


EARLY 

THERMAL FITS TO 

CAS A AND 

TYCHO 



Cas A 


Tycho 


Ref 

Tlow 

Thigh 

Tlow 

T high 

Gorenstein et al. 1970 

— 

1.25± .34 

M •> — M 

2. 2*1.0 

*Brisken 1973 

. 8 J± . 06 

- - - - 

.45*. 05 


Hill et al. 1975 

— 

2.3 

— 

1.9 

Charles et al. 1975 

.73*. 20 

2.6 ±1*0 

— 

— 

Pravdo et al. 1976 

< *7 

3.9 ± .5 


— _ 

Davison et al. 1976 

1.08*. 03 

5.2 * .6 

.55*. 04 

3.5* .3 


^refined analysis of data of Serlemitsos et al. (1973) 


B. THE EINSTEIN SSS DATA 

The newer Einstein data summarized here are from the Solid State 
Spectrometer (SSS), a cryogenically-cooled Si (Li) device with a FWHM 
energy resolution of ~ 160 eV over its effective 0.5-4. 5 keV energy 
range (see Holt et al. 1979 for experiment details). The spectra I 
shall discuss have already been prel iminarily reported in the 
literature. The new information presented here will be of two kinds: a 
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refined set of experimental spectra with all sources of background over 
the full 0.5-4. 5 keV band now properly accommodated, and some remarks 
about the progress which has been made in the physical modelling of SNR 
to reproduce these spectra. 

A Crab spectrum is presented in Fiqure 1; it represents one 
of ~ 100 such spectra obtained over the lifetime of the experiment, all 
of which were used for calibration purposes to ensure uniformity in the 
interpretation of detector response with time. The input spectrum is 
presumed to be a featureless power law of index a - 1.1 observed through 
an interstellar medium of cosmic abundance with column density 
N h * 3 x 10 21 atoms cm-2. Note that the experiment introduces *,n 
apparent kink in the spectrum at ~ 0.8 keV {arising from the nickel with 
which the telescope is coated), but that any other features 
(particularly ones assocaited with silicon, of which the detector is 
made, and which might be expected near 1.8 keV) are much smaller than 
the statistical precision of one of the displayed pulse-height channels 
in a typical measurement. The Crab spectrum is illustrative, therefore, 
of the fact that any potential line features in the spectra are not 
systematically induced in the detector. 

The line features which are observed in young SNR spectra are 

indicated in Figure 2, where are plotted the pulse heioht data from an 

exposure to the Tycho (SN 1572) remnant, toqether with a comparison 
spectrum for illustrative purposes. The comparison spectrum is that 
expected from a T = 0.5 keV solar-abundance plasma in thermal 

equilibrium, viewed through an interstellar column density of 2 x 10 21 

atoms cm-2; this comparison histogram is not convolved throuqh the 
detector response, but is smeared to a FWHM resolution of ~ 50 eV with 
exactly the same bin widths as the experimental data. The shaded 
portion of the comparison spectra is the contribution from Fe-L emission 
from the isothermal plasfta, and the darkened portion the contribution 
from Si-K emission; the four blackened bumps represent 
Si +12 Ka, Si +13 « aj si +12 and Si +12 Ky in order of increasing 
energy. From the obvious line identifications, therefore, it is clear 
that the observed lines are characteristic of a relatively low 
temperature (e.g. T = 0.5 keV) plasma, in the sense that the helium-like 
lines of the Z > 10 constituents are evident while the hydrogen-like 
lines are not, hut the equivalent continuum widths of the observed lines 
are considerably in excess of the equilibrium values for solar 
abundances. Most of the remainder of this paper will be devoted to a 
discussion of the model! inq of such spectra. 

C. SPECTRAL MODELLING 

Figure 1 displays the SSS spectrum from Cas A. As originally 
reported by Becker et al. (1979), the data cannot be fitted with a 
single equilibrium temperature with any abundance distribution; instead, 
it can be reasonably well -fitted with two equilibrium components. The 
lines are consistent with an equilibrium plasma at T - 0.65 keV (albeit 
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Figure 1. Experimental SSS spectra for the Crab nebula, CasA and 
SN1Q06. The solid histogram through the SM1006 data points represents a 
power law of index a = 1.2 observed through a column density of 9 x 1020 
atoms cm-2. 

with overabundances in the line-emitting components), but a higher 
temperature component is required to replicate the harder continuum 
observed and to smoothly connect to the hioher enerqy measurements; 
because the bremsstrahlung conti nua from any temperatures £ 3 keV are 
indistinguishable in the SSS 0.5-4. 5 keV bandpass, a fixed 4 keV was 
chosen to approximately connect to the higher energy measurements. 

The data presented in Figure 1 are somewhat more refined than 
originally presented in Becker et al. (1979), but the important 
source-fitting parameters are not substantially changed. In particular, 
the best-fit low temperature component (for a 2-temperature equilibrium 
fit) is still T « .65 keV, and the Si abundance is ~ 1.5 times solar, as 
determined from the simultaneous fitting of the lines and the 
continuum. Relative to this inferred Si abundance, the inferred 
abundances of other important constituents relative to solar proportions 
are S/Si{~2), Ar/Si{~4), Ca/Si(~2) and Mg/Si (~.l) . The only significant 
difference between the results of the present modelling and that 
originally reported is in the Fe abundance, as the model extension to 
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ENERGY (keV) 

Figure 2. Experimental SSS spectrum of the Tycho SNR, compared with the 
expected spectrum from a solar-abundance equilibrium plasma at T « 0*5 
keV as it would be observed with a detection system which has an 
energy-independent response (see text for additional details). 

lower energies changes what was previously a measured Fe/SI ~ .5, 
relative to solar, to an upper limit. The lack of a pronounced increase 
with Z in the relative overabundances of Si, S, Ar and Ca suggests that 
the standard qualifications regarding non-equilibrium ionization may 
obscure the important qualitative conclusions that may be drawn from the 
spectrum. If, as suggested by Fabian et al. (1980) from an 
interpretation of the X-ray morphology, the remnant is in the free 
expansion phase and the low energy X-ray emission is dominated by the 
ejecta, the dumpiness of the emission suggests a situation closer to 
equilibrium (at least within the clumps) than would a smoother 
distribution, as the Rayleigh-Taylor-induced density enhancements should 
equilibrate more rapidly than would the unclumped ejecta (see Chevalier 
1982 for a discussion of the unclumped' density distribution which might 
be expected from a system like Cas A with pre-SN mass loss). The modest 
abundance enhancement of Si-qroup material implied by the equilibrium 
simplification may, therefore, better characterize the spectrum than 
would a detailed (non-unique) calculation, at least until such time as 
measurements with better spectral and spatial resolution are available. 

In contrast to Cas A, the three Type I remnants Tycho, Kepler and 
SN1006 display a clear necessity for non-equilibrium effects. The 
spectra of Tycho and Kepler bear marked similarities to each other (see 
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Figure 3); utilizing the same two- temperature equilibrium model as a 
starting point, both have low temperatures of T « .5 keV. Tycho 
requires a Si abundance which is approximately an order of magnitude in 
excess of solar (Kepler requires about half the Tycho Si abundance), and 
for both remnants the allowable relative abundances of S, Ar and Ca 
increase with atomic number (in the newer equilibrium fits, these 
implied abundances are even higher than those in Becker et al. 1980a and 
Becker et al. 1980b), in accordance with the expectation for a 
non-equilibrium plasma which is still ionizing. The application of a 
Sedov approximation (Shull 1983; Hamilton, Sarazin and Chevalier 1983) 



Figure 3 - Experimental SSS spectra of the Tycho and Kepler SNRs. The 
solid histograms, arbitrarily placed above and below the Kepler 
spectrum, represent different fits to Kepler of comparable statistical 
adequacy. The upper fit is a two-temperature equilibrium model for 
which the abundances are free parameters (see text), while the lower fit 
is a Sedov model with similarly unconstrained abundances; many such 
models (with differing q, electron fitted with comparable statistical 
adequacy, but all require overabundances relative to solar. 
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relieves this latter problem, but its applicability may be 
problematic, A variety of Sedov models with and without clumping (and 
with and without the assumption of rapid equilibration of the electrons 
with the shock) all yield pronounced overabundances of Si, comparable to 
those found in the equilibrium fits. The relative abundance ratios for 
the higher-Z constituents (i.e. those which are identifyable by their 
K-emission) are no longer unreasonably high, however; S/SI is still ~ 2 
(as it appears to be for all SNR fits), but Ar/Si and Ca/Si need not be 
much higher than unity (i.e. solar). The utilization of a Sedov model 
(properly constrained to the observable luminosity at the best-estimated 
distance, with reasonable values of n 88 E 0 n 0 2 ) is questionable a priori 
if the interpretation of the X-ray morphology given by Seward, 

Gorenstein and Tucker (1983) is correct, i.e. if the large fraction of 
the X-ray luminosity arises from clumped reverse-shocked ejecta rather 
than from the blast wave. Interestingly, Hamilton (this volume) has 
demonstrated that the emission from the ejecta can be modelled with a 
Sedov approximation (where the parameter n is a bit more obscure), 
allowing fits which better match the data in detail (Sarazin, Hamilton 
and Chevalier, this volume). 


A continuing feature of all the fits is the lack of marked 
overabundances of Fe, as might be inferred from its L-emission (see 
Figure 2) in accordance with current models of Type I SN which release a 
large fraction of a solar mass of Fe-group material. The primary 
difference between the Tycho and Kepler spectra of Figure 3 is in the Fe 
L-emission; Fe is virtually undetectable from Tycho, but is required at 
about the level Fe/Si ~ 1 (relative to solar) in the various equilibrium 
and non-equilibrium Kepler fits. Even in Kepler, this abundance is at 
least an order of magnitude below that expected if the Fe-group material 
is well -mixed into the X-ray-emitting plasma. Clearly, the retention of 
the Fe-rich model requires the not-unreasonable assumption that the 
heavier ejecta is interior to the Si -group ejecta, and has not yet been 
si gni f icantly reverse-shocked. 

The X-ray spectrum of SN1006 is even more puzzling. Unlike the 
younger Type I remnants discussed above, there is no statistically 
significant evidence for K-line emission. In fact, the X-ray spectrum 
can be fit with a Crab-like power law (a ~ 1) from 0.5 keV out to at 
least 20 keV (Becker et al. 1980c). It is worth noting that the higher 
energy X-ray spectra of other SNR can be relatively well-fit by power 
law approximations, but with spectral indices which are much larger 
(i.e. the younger and more Initially energetic SNR Cas A can be fit 
with a ~ 3, while Tycho can be fit with a ~ 4; see Figures 1 and 2 of 
Pravdo and Smith 1979). Fabian, Stewart and Brinkmann (1982) have 
suggested that the appearance of a featureless power law may arise from 
the superposition of emission from differing thermal shells. If the 
intermediate-Z metals are substantially locked in grains, and if the 
reverse-shock has long-si nee passed through the Fe-rich ejecta (so that 
the remnant interior has cooled to ~ 11)6 K), the absence of the 
prominent lines observed in the younger remnants may be reconciled. It 
should be noted, however, that the spectral hardness of the high energy 
continuum may be more easily explained with non-thermal processes, e.g. 


20 


ORIGINAL PAGE W 


synchrotron radiation from electrons accelerated near the blast wave 
(Reynolds and Chevalier 1981) or possibly even from the naive picture of 
a hidden pulsar which was so popular ten years ago. 


D. FINAL REMARKS 


The combination of higher-energy (> 5 keV) spectral data from 
pre-Einstein Investigations (see Pravdo, this volume) with morphological 
and spectral data from all the Einstein instrumentation provides a 
powerful tool for studying SNRs. This evidence is not totally 
conclusive in all particulars, however. The modelling of the SSS data 
requires a complete time-dependent coupled solution of the hydrodynamics 
and ionization, including whatever density inhomogeneities may exist in 
each specific remnant, since the experiment aperture is 6 min in 
diameter. Better spectral resolution than its 160 eV FWHM would be 
useful (insofar as K-emission contributions from differing multi pi et 
components would then be resolvable), but the present resolution is 
sufficient to distinguish H-like and He-like ionization states. Much 
more important would be the application of the available spectral 
resolution to much more localized regions of the SNR, so that 
differential information would not have to be extracted from the 
attempted separation into components of the whole-remnant spectrum. 

The available procedure is necessarily non-unique, and cannot 
unambiguously distinguish whether the electrons in the recently shocked 
plasma can quickly equilibrate with the heavier constituents (as it 
least some fraction must, from > 10 keV measurements; Pravdo and Smith 
1979) or whether they require a Coulomb equilibration timescale 
(typically longer than the remnant age) to do so; the two cases yield 
very different bremsstrahlung conti nua and, hence, can result in very 
different blast/ejecta proportions and abundances. Unfortunately, fits 
to the whole remnant data with either assumption for the electron 
equilibration timescale can yield results of comparable statistical 
significance. With no independent reason for fixing such otherwise 
indeterminate parameters, therefore, model uniqueness is not assured. 
Until such time as more independent information exists (e.g. from the 
application of like quality spectral data with much finer spatial 
resolution, so that some separability of otherwise degenerate fitting 
parameters can be obtained), a satisfactory fit to the whole-remnant 
spectral data cannot guarantee a unique detailed explanation. 


My thanks to the many members and associates of the GSFC X-ray 
group, past and present, who have contributed to our ongoing program of 
research. Andrew Szymkowiak, in particular, is responsible for the 
large fraction of the previously unpublished work on which much of this 
discussion is based. 
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Kirshner: The principal optical result on Cas A is that oxygen is 

greatly enhanced. The SSS spectra only give the ratio of Si, S and the 
rest to oxygen. Are the X-ray data consistent with the extraordinary 
abundances seen in the optical filaments? 

Holt: Because the energy resolution of the SSS does not allow us to 

completely separate line and continuum contributions to the experimental 
spectrum, and because the bremsstrahlung continuum is dominatted by 
oxygen (if present in at least solar proportions to hydrogen), it is not 
possible to define a totally unambiguous reference for the 
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"abundances' 1 — even before addresslnq the problem of how appropriate an 
equilibrium model should be for the whole-remnant specrum of Cas A. If 
we assume that the portion of the spectrum below ~ 3 keV Is dominated by 
the ejecta in the fast-movinq knots, the implied Si/0 ratio is well 
within a factor of two of the values given in Table 5 of Chevalier and 
Kirschner (1970). 

Tuffs: Is a shock velocity of, say, 2800 km s _1 for Cas A consistent 
wTfFHthe observed line ratios? 

Holt: The shock temperature defined by the shock velocity (for a Sedov 
solution) is > 1QB K, or much too hi qh to match the observed line 
ratios. The experimental preponderance of helium-like to hydrogen-like 
lines forces the equilibrium fits to the spectra of Tycho and Kepler, in 
addition to that of Cas A, to contain a component at < 10? K which 
dominates below ■« 3 keV, At higher energies, however, there is 
continuum evidence for plasma approachinq equilibrium with the blast 
wave (Pravdo and Smith 1979). 

Chevalier : With regard to the non-uniqueness of the X-ray spectral 
T-'Tts 71<T a! 1 the models have the proper ages, radii and X-ray 
luminosities? 

Holt: The Sedov models we have attempted to fit to Tycho and Kepler, 
“such” as that displayed in Figure 3 for Kepler, are those for which the 
ages are correct and the radii /luminosities are consistent with 
reasonable distance estimates. We don't take seriously any fortuitous 
spectral fits with unreasonable model implications. 

Fedorenko: Are the power-law X-ray spectra consistent with the 

equi librium particle distribution functions? What is the mechanism of 

radiation? 

Holt : I presume that you are referring to the spectrum of SN1Q06, for 

which we find that a power-law is an 'adequate fit to our data. It is 
certainly possible to synthesize power-law fits from combinations of 
equilibrium particle distribution functions (see the early attempt by 
Sartori and Morrison 1967 to so reproduce the spectrum of the Crab), but 
the consistency of the same power law over two orders of magnitude in 
enerqy is suqgestive of a synchrotron origin for the X-ray emission. 

Fabian; Since the electron-ion equilibration time is often greater than 
the age of the remnant, how do you know what high enerqy spectral form 
should be fit to the continuum data (e.g. in SN1 006 ) ? 

Holt: In the case of Sedov models, we utilize computer-generated 

differential spectra instead of limiting analytic forms, which are 
dependent upon the usual model parameters (including time) as well as 
upon physical assumptions such whether electron-ion equilibration can be 
more rapid than specified by Coulomb interactions alone. Althouqh a 
power-law appearance is certainly possible for the superposition of 
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non-equilibrium as well as equilibrium thermal spectra, the extension of 
the form to at least 20 keV for SN1Q06 argues against it because such a 
hard spectrum (a ~ 1) suggests a blast wave with much higher velocity 
than that of the younger and more energetic Gas A, for example 
(with a <v 3); the appearance of any photons at all at 20 keV requires a 
current blast wave velocity of 4000 km s-1, but the maintai nance of the 
hard spectral form requires that the blast wave velocity is considerably 
larger than that. The assumption of a slow equilibration time only 
exacerbates the problem, as it is then even more difficult to produce 
X-rays with energies comparable to the value for equilibrium with the 
blast wave and, hence, more difficult to produce both a hard spectrum 
and its extension out to 20 keV, 
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ABSTRACT 

We present soft X-ray photomosaic Images of two supernova 
remnants, Puppls A and IC 443, constructed from a series of exposures by 
the Einstein Imaging instruments. The complex morphologies displayed In 
these Images reflect the Interaction between "middle-aged" Bupernova 
remnants and various components of the Interstellar medium. Surface 
brightness variations across Puppls A suggest that Inhomogeneitiea on 
scales from 0.2 to 30 pc are present in the Interstellar medium, while 
the structure of IC 443 is apparently dominated by the Interaction 
between the remnant and a giant molecular cloud. 


In supernova remnants that have evolved well Into their adiabatic 
expansion phase (t ** 10^ yr), the bulk of the soft X-ray emission arises 
from recently shock-heated interstellar matter. The soft X-ray 
morphology of such "middle-aged" remnants will be influenced by any 
large or small scale inhomogeneitles In the interstellar medium (ISM) 
they encounter. Consequently, these supernova remnants serve as 
excellent remote probes of the structure of the ISM. 

We present and briefly discuss here high-resolutlon soft X-ray 
images of two such "middle-aged" remnants, Puppls A and IC 443. Despite 
their similarity in age (~ 10^ yr), distance (~ 2 kpc) and diameter 
(25-30 pc), these two remnants display markedly different morphologies, 
presumably due largely to their different environments. A more 
extensive discussion of Puppls A may be found in Petre et al. (1982a); 
the IC 443 Imaging results, along with spectral studies using the 
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Figure 1 ; A high-resolution 0.1-4 keV photomosaic map of Puppls A. 

Einstein Solid State Spectrometer (SSS) and the HEAO 1 A-2 MED, will, 
appear in Petre et al. (1982b). 

Puppls A 

Figure 1 shows a high-resolution X-ray image of Puppls A in the 
energy range 0. 1-4.0 keV. The map is an exposure-corrected photomosaic 
of 11 Einstein High Resolution Imager (HRI) exposures, binned in 8" 
pixels. The image reveals an ISM with structure on many scales. First, 
there is a general decrease in X-ray surface brightness, by a factor of 
at least 20, from northeast to southwest across the remnant, 
perpendicular to the galactic plane. This effect is probably due to a 
local density gradient in the ISM of at least a factor of 4 over a scale 
of ~ 30 pc. Second, the surface brightness variations in the interior 
and along the shell suggest a preshock ISM containing many small density 
variations (factors of 2), about an average of " 1 cm , over scales 
from 0.2 to 10 pc. Despite the surface brightness variations along the 
shell, the pronounced limb -brightening profiles are for the most part 
consistent with that expected for an adiabatlcally expanding blast wave, 
suggesting that the local dumpiness of the interstellar matter docs not 
prevent the adiabatic model from being locally valid. Finally, in 
addition to the filamentary structure representing small density 
perturbations, Puppls A contains two bright knots of emission along the 
shell, one in the east and one in the northwest. These knots are 
probably shocked Interstellar clouds with preshock density of 10-20 
and diameter 1-2 pc. 
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Figure 2 deplete e 0.2-3. 1 keV map of IC 443, overlayed on a 
Palomar Sky Survey red plate. The map la an expoaure-corrected , Wiener- 
filtered photomoaalc of three Imaging Proportional Counter (IPC) lmagea. 



Figure 2 (left): A 0.2-3. 1 keV Image of IC 443, auperpoaed on a Palomar 
Sky Survey red plate. Thia map la a photomoaalc of three IPC 
lmagea. F igure _3 (right): An HRI photomoaalc of IC 443, auperpoaed 
on the Palomar red plate. 

binned In 32" plxela. The contoura scale linearly with surface 
brightness. Although X-ray emission la detected from the entire SNR, 
the remnant does not appear to be llmb-brlghtened. The average surface 
brightness in the vicinity of the bright northeastern optical filaments 
la a factor of * 4 higher than elsewhere. Thia network of filaments has 
apparently been created by the collision of the shock front with a giant 
molecular cloud, which la observed In CO to extend well beyond IC 443 to 
the northeast (Scovllle et^ a_l . 1977; Cornett, Chin and Knapp 1977). 

A sharper view of the northern portion of IC 443 la shown in 
Figure 3. This map is an exposure-corrected, Wlener-flltered 
photomosnlc of three HRI Images, with 16" binning, whose coverage of IC 
443 Is represented approximate ly by the area over which the grid 
appears. The contour levels w- re chosen to reveal the highest X-ray 
surface brightness features; X- «ys were actually detected by the HRI 
from everywhere In the complex <f optical filaments, with the notable 
exception of the finger of opt al extinction extending northeast ward 
into the filaments. 

The HRI map emphasises he dominant X-ray feature, a * 15' 

diameter region whirl partiall overlaps the optical filaments. The 
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pressure of the X-ray emitting gas in this region, as inferred from 
spectral measurements by the Einstein SSS, is ~ 3 x 10“” dyne cm" 2 , 
placing it in approximate equilibrium with the cooler, denser, 
optically-emitting gas (Feseu and Klrshner 1980). Coupled with the 
general overlap between the regions of high X-ray surface brightness and 
the bright northern optical filaments, this suggests that the X-ray 
emitting gas may be interspersed among the sites of unstable radiative 
cooling which appear as optical filaments, and that both have arisen 
from the collision of the shock front with the molecular cloud. 

The molecular cloud may affect the observed low-energy structure 
of 1C 443 in a second way. As is visible in Figures 2 and 3, a lane of 
high optical extinction bisects IC 443 from northwest to southeast, 
extending a finger into the bright optical filaments. This lane 
demarcates another portion of the molecular cloud. The detection of 
shocked [H I] and a va~iety of molecular species with typical velocities 
around -30 km s” 1 suggests that the front of the IC443 shell is 
colliding with the molecular cloud (see, e.g., Giovanelll and Haynes 
1979; DeNoyer 1979). Although no general correlation between the X-ray 
surface brightness, as measured by the XPC, and the column density of 
the cloud, as Inferred from the CO brightness temperature, is apparent, 
a correlation between column density and X-ray spectral hardness ratio 

nrnml a* 


Figure 4,: An IPC spectral hardness 

ratio map of XC 443, superposed 
on CO brightness temperature 
contour® (from Cornett, Chin 
and Knapp 1977). Spectral 
hardness is defined as the 
ratio between 0.2-1. 3 kev 
counts and 1.3-3. 1 keV counts, 
with symbols of increasing 
prominence representing ratios 
of 0.6, 0.8, 1.0 and 1.2* 

Circle is HPBW of CO map. 

6" 15*0' 6*MV 6 h |3"0' GW 

(0.2-1, 3 kev/1. 3-3.1 kev) does exist. As seen in Figure 4, X-rays from 
regions coincident with the cloud are harder (average hardness ratio of 
~ 0.8) than elsewhere in the remnant (average hardness ratio of ~ 0.5). 
To account for this increase of spectral hardness by absorption, a 
column density through the cloud of 2-4 x 10 21 cm” 2 is required (in 
addition to the column density to IC 443) if the characteristic 
temperature of IC 443 is uniform and within the range of 0.5-1. 5 kev. 
This additional column density may be compared to the ^inf erred average 
column density of H2 through the cloud of 2 x 10 21 cm” (i.e., 4 x 10 2 
cm” 2 of [H I] - Cornett, Chin and Knapp 1977). It is thus probable that 
the cloud absorbs X-rays emitted by the interior of IC 443. 
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JL* Introduction 

In this presentation we show how the study of the isotropy of the 
X-ray sky contributes to our understanding of the structure of the 
universe at moderate redshifts (l^ z<<z re combination^ * Actually, the 
anisotropy of the sky flux provides the information, much as the 
microwave sky anisotropy does for earlier epochs. [See reports in this 
volume.] Though we are currently unable to make measurements with the 
precision and small solid angles typically achieved in the microwave, 
comparatively crude limits from the X-ray fluctuations place limits on 
the largest scale structure of the universe. We first outline the 
measurements of the X-ray sky and its anisotropies made with the HEAP 1 
A-2 experiment. Detailed presentations are found elsewhere [Shafer 
1982; Marshall et al. 1980; Piccinotti et aJL. 1982; Iwan et_ al . 1982; 
Shafer et_ al. in prep.]. We then show how the anisotropies place limits 
on the origin of the X-ray sky and on any large scale structure of the 
universe, following the example of previous analyses which used earlier 
anisotropy estimates [see _e.&. Fabian and Rees 1978; Rees 1980; Fabian 
1981] . 

2 . The X-ray Sky 

In Figure 1 we present the extragalactic sky spectrum. Several 
properties of the X-ray portion of the spectrum are noteworthy: 

(1) It is bright. A spectrum with slope -1 in Figure 1 has equal 
energy per decade; thus the energy density in 3-100 keV X-rays is second 
only to the density of the microwave region. 

(2) It Is easily detected and nearly isotropic. The region from about 
3 keV to ~1 MeV Is the only well studied portion of the sky spectrum 
other than the microwave background that is not dominated by a strong 
galactic component. 

(3) It has a well determined, if not well understood, spectrum. In 
terms of accuracy and bandwidth the measurement of the X-ray spectrum 
surpasses even determinations of the microwave spectrum [de Zotti 
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Figure 1. The Cosmic 
Extragalactic Spectrum. The 
light dashed line extends the 
radio discrete sources to 
higher energies assuming a 0.7 
energy spectral index. [Radio 
and microwave: Longalr 1978; 
Optical: Dube, Wickes and 
Wilkinson 1979; UV: Paresce, 
McKee and Bowyer 1980; X-ray: 
Marshall .et al. 1980, 
Rothschild jet al. 1983; y-ray: 
Flchtel & Tromblca 1981] . 


1982a] . The spectrum from 3 to 100 keV is well represented by a thin 
thermal bremsstrahlung model with a temperature of 40^5 keV [Marshall elt 
al . 1930; Rothschild et al. 1983]. Mo population of sources with a 
single index power law spectra provides the right shape. A suitably 
constrained model of an evolving population with a sharp break in their 
spectral index is consistent with the 3-50 keV data [de Zotti et 
al. 1982b]. 

A variety of models have been presented to account for parts of 
the sky flux, but no coherent picture exists that satisfies all the 
observations. Particular models may be consistent with the data but all 
leave open questions or await particular observational confirmation. 
Known and possible fractions of the X-ray sky flux are: 

(1) A galactic component. Because of the observed Isotropy of the 
sky, any flux associated with the galactic disk can contribute no more 
than ^2-10% of the high galactic latitude flux, depending on the 
latitude of the observation [see e_.£. Iwan ejt_ a_l. 1982] . In addition, 
known populations of galactic sources do not generally have the same 
spectrum as the sky. 

(2) Well observed extragalactic sources. Based on the HEAP 1 A-2 all- 
sky survey the extragalactic sources resolved at high flux are 
predominantly either clusters of galaxies or active galactic nuclei 
(Seyferts, N galaxies, etc.) [Piccinotti et. al. 1982]. The derived 
local luminosity function of these two populations can be used to 
estimate contributions to the total sky flux of about 4% and 18% 
respectively, assuming no significant evolution of the luminosity 
function. In addition, the low temperature thermal spectra of cluster 
sources [Mushotzky jet al . 1978] and the power law spectra with a photon 
index near 1.7 typical of active galaxies do not correspond to the sky 
spectrum [Mushotzky e£ al . 1980; Rothschild jet al. 1983]. In fact, no 
source or population has been observed to have the proper spectrum to 
provide the bulk of the sky emission. 

(3) An evolving population of sources. This could be the evolution of 
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the above known populations or the Introduction of a new population of 
sources. QSOa, which are undergoing apparent evolution at other 
wavelengths, have been shown to be strong emitters in the softer 0.5 to 
3,5 koV band covered by the Mnstyoljy observatory Csee &.&» Xamoranl tvfc 
nl» 19811 « Unfortunately, there are too few high quality broad band X- 
ray spectra of QSOa to make an unambiguous estimate of their 
contribution at the higher energies typical of the bulk of the X-ray sky 
flu*. Aval [1078] pointed out that active galaxies may make up the 
total sky llux if they undergo only moderate evolution, in comparison to 
the amount ol evolution suggested for QSOs in the optical. However, the 
evolution must involve the spectral form of the objects as well as their 
luminosity function [tyjfi. belter and lkildt 1982]. Suggestions for new 
populations of X-ray sources have included hot gas associated with the 
initial generation of stars of young galaxies (bookbinder et, al. 1980] 
and primordial black boles (Garr 1980], Though the proposed spectra are 
in accordance with the X-ray sky spectrum, there have been no 
identifications of these new objects with observed X-ray sources, 

(3) A totally diffuse component, such as a hot intergalactic medium. 
This model baa the correct spectral form, but there are possible 
difficult ten providing the energy to heat the medium (field and Perronod 
1977,* fab lan 1981], 

Though the exact origin of the sky flux in still an open question, 
we have a better understanding of the principal sources of the observed 
anisotropy* We classify the variations in the X-ray sky intensity as 
large angular scale anisotropies or as fluctuations (small scale 
variations') , 

The galaxy does not dominate the 2-10 keV sky flux, but the 
dominant large scale varia tion is associated with the galactic disk. An 
early model of this variation was the cosecant |b| law of an infinite 
plane of emission (Warwick, Pyo and Fabian 1980]. Other studies have 
noted a longitudinal component associated with the galaxy [Protheroe, 
Wo 1 feudal e and Wdoweayk 1980; Xwan et, aU 1982]. An expected large 
scale anisotropy, of smaller magnitude, is a cosine or dipole 
anisotropy. Such a signal is expected for the. same reason as the dipole 
variation seen in the microwave sky, i..,o. motion of the observer with 
respect to the rest frame of the emission, the Compton-Got ting effect. 

The dominant contribution to the small scale fluctuations is a 
continuation of known source populations to lower flux levels where 
sources are no longer individually detectable. The. also and shape of 
the frequency distribution of the fluctuations is a function of the 
number of sources versus flux relation, N(S), at those fluxes. Wo can 
explain all the fluctuations in terms of the known populations, without 
evolution, and place an upper bound on the size of any other small senle 
varlations, hereafter referred to as the excess variance. This bound 
constrains all other sources of anisotropy. Possible origins of 
additional variation would bo an evolving or new population of sources, 
or a clumping of the sources that make up the background. At the 
largest angular scales source o lumping indicates large scale structure 
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In the universe, such as a global perturbation In the density Sp* An 
upper bound oq_the excess variance limits the allowed strength of such 
structure, ®P/P, with no assumption about the origins of the X-ray sky 
flux other than presuming that variations in the X-ray volume emlosivlty 
are proportional to <Sp. (For general reviews of the X-ray sky see c*ji. 
boldt [1981], Fabian [1981].) 

3. The Data 


Our results are based on measurements taken with a xenon 
proportional counter, one module of the A- 2 experiment on the HEA P; _l t 
satellite [Rothschild ejt aJL. 1979] , taken during an all-sky survey* For 
the X-ray sky spectrum, 90% of the counts in this detector originate ^,n 
the 2.5-13.3 keV band. We measure flux, S, in units ^of counts s" 1 cm"". 
For typical extragalactic spectra 1 count a" 1 cm"" is equivalent to 

I. 35xl0" 8 ergs a"* cm"** (2-10 koV). The all-sky flux, S as , is 58 counts 
s" 1 cm . The measured count rata depends on detector area, collimator 
solid angle^ and integration time. For our measurements the mean sky 
intensity, l s j. , is 17.06 counts exp" 1 (on* 4 exposure is 1,28 a). The 
mean count rate of the internal, non-X-ray, background is 3.5 counts 
exp” 1 . The average uncertainty due to counting statistics was 0.23 
counts exp . The angular size of the measurements is fairly large, 
over 100 square degrees, but much of this urea contributes little to the 
total count rate* 90% of the total comes from a rectangle of 

II. 2°x4.4 0 , covering 49 square degrees. The central area of ~26 square 
degrees contributes 71% of the sky intensity but 90% of any excess 
variance in the intensity. 

The fluctuation data were restricted to high galactic latitudes, 
jb|>20°, and free from contamination by X-ray sources in the Magellanic 
Clouds and bright high latitude galactic sources, When looking for 
large scale structure, we included data down to latitudes of 10°, 
excluding all contamination from any resolved source cataloged in the 
complete all-sky sample of Plccinotti et al* [1982] • 

The A-2 detectors had several unique features for the. continuous 
monitoring of internal background and determining the X-ray sky flux. 
The performance of the detectors, as monitored by repeat scans of the 
same area of the sky six months apart, was very stable. The internal 
background was also very stable. After selection of data to avoid noisy 
periods, the variation in the background was roughly 0.05 counts exp" 1 , 
corresponding to a sigma 1.3% of the non-X-ray count rate and only 0.25% 
of the total intensity. 

A* barge scale variations : Galactic and Dipole anisotropies 

Using a similar set of HEAP 1 A-2 data, Iwan .et ol- [1982] showed 
that there was a variation associated with galactic longitude in 
addition to the latitude variation. Following that paper we fit: the 
galactic component with a disk of finite radius and an exponential scale 
height. These parameters are strongly correlated and their upper bounds 
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poorly determined. The best fit radius was 1.8 Rg C , where Rg C Is the 
distance from the sun to the galactic center, roughly 10 kpc. The best 
fit value for the scale height ms 0.4 R gc » with the 90% lower limit of 
0.lR gc , corresponding to 1 kpc, a scale typically larger than most 
galactic X-ray source populations. A discussion of this model and its 
implications is given in Iwan et al. [1982]. 


After removal of the best fit calactic model, we fit a dipole 
model, SI « I GG cos 0, where 0 is the angle between the observation and 
the signal maximum. The addition of this new model to the fit produced 
a drop in x 2 significant at the 95% level. The strength of the signal, 
I CG , is 0.0910.03 counts exp" 1 , about 0.5% of the sky intensity. The 
best fit direction in galactic, coordinates, (£,b), is (282°, +30°) . 
However the 90% confidence region for the direction is very large, 
covering about one eighth of the total sky. A result of similar 


direction, magnitude, precision, an 
Wolfendale and Wdoczyk [1980] using 


Figure 2. Position of Dipole 
Maximum. Contours show 70% 
and 90% confidence regions. 
The center + marks the best 
fit position. Also shown are 
the one-sigma error bars for 
measurements of the dipole 
maximum in the microwave. 
[BCW; Boughn, Cheng & 
Wilkinson 1981; FGMN: Fabbri, 
Guidi, Melchiorrl and Natale 
1980; GSs Gorenstein and Smoot 
1981.] 


1 confidence was found by Protheroe, 
UHURU data. 



One possible interpretation of this statistically marginal result 
is in terms of the Compton-Getting effect, where the size of the dipole 
signal is related to the observer's velocity by 

I CG * I (2+r) v/c. [1] 

P is the photon index of the sky flux, ~l-.4 for the band we are 
interested in. The derived value for the velocity, 475+165 km s" , and 
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the direction are consistent with observations of the dipole signal in 
the microwave sky. The microwave directions in the literature are shown 
in Figure 2 along with the X-ray confidence regions. A synthesis 
performed by Wilkinson at this symposium of the different experimental 
results gave a best fit direction of (265°»+50°) and a velocity of 
372±25 km s~ . The origins of the Compton-Getting velocity may be 
responsible for an additional component of the observed X-ray dipole 
signal. If the velocity observed in the microwave is the integral of 
the acceleration caused by a large scale overdensity ("lump") at the 
same redshlfts range at which the X-ray sky emission originates, then 
the lump should produce an excess in the X-ray emlsslvity, producing an 
enhancement in the direction of the lump in addition to the Compton- 
Getting velocity signal. Comparisons of the X-ray and microwave large- 
scale anisotropies help to decouple the two effects of such a lump, 
providing constraints on the lump's properties. [see js.£. Warwick, Pye 
and Fabian 1980; Fabian 1981]. 

The total dipole signal may be due to variations other than the 
Compton-Getting effect. The form of the underlying structure may not be 
a pure dipole. For instance, there may be second order anisotropies 
associated with the galaxy but not included in the finite disk model. 
Also, it is intriguing that the contour for the direction of the dipole 
maximum includes a large fraction of the local supercluster, which could 
provide a large scale enhancement to the sky flux. The best fit value 
of Iqq corresponds to a maximum surface brightness of 0.02 counts s“ 1 
cm”^ sr“ . The volume emissivlty and total luminosity of the local 
supercluster required to dominate the "dipole" signal are dependent on 
geometry. Assuming a disk of emission centered on the dipole maximum 
with a radius of 4 Mpc, an estimate of the total luminosity is 2x10^ 
erg s - ” 1 (H 0 = 50). Previous attempts to correlate X-ray surface 
brightness with the local supercluster have yielded only upper limits 
larger than the above estimates. [see £.£. Schwartz 1980]. To test if 
the local supercluster is in part responsible for the fit dipole signal 
requires direct testing of models for the supercluster, a project now in 
progress. The tidal "12 hour" signal reported by Warwick, Pye and 
Fabian [1980] in the Ariel V data was not observed in the A-2 data. 

5.. Small Scale Anisotropies : The Fluctuations 

Point sources have an impact on measurements of the sky flux even 
if the sources are too numerous to be individually resolved. The actual 
number of sources, and hence their total intensity, varies from one part 
of the sky to another. The process of extracting information about the 
sources from the size and shape of the intensity distributions was 
pioneered by radio astronomers [Scheuer 1957; Condon 1974; see 
j 2 .&. Condon and Dressel 1978]. Given a model for the differential 
number of sources as a function of flux, N(S) dS, the distribution of 
intensities, Pj_y(I) dl, can Predicted, assuming the distribution of 
sources is completely random and unclumped. Standard statistical tools 
are used to evaluate the N(S) models by comparing the predicted 
distributions to observations, extending our knowledge of the X-ray 
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source counts beyond what was directly accessible from resolved sources 
[Fabian 1975; Schwartz 1976; Pye and Warwick 1976]. 


Figure 3. Number of Sources 
Versus Flux Derived from the 
Fluctuations. The plotted 
quantity is the ratio of the 
number to that expected for a 
fiducial Euclidean model, 
where N * S”” 5 ' • See text for 
details. 



In Figure 3 we present the results of such model fitting using the 
HEAP 1, A-2 data. We restricted our model to a single power law form, 

N(S) dS * 4ir K S~ Y dS , [2] 

with a sharp cutoff imposed at the low flux where the intensity 

contributed by the sources equals the total sky flux: 

/ dS N(S) S » S us . [3] 

A population of sources distributed randomly through Euclidean space 
will follow a power law model with Y - 5/2. We compare an N(S) relation 
to the Euclidean form by defining a function K(S) 

N(S) 5 4ir K(S) S~ 5/2 . [4] 

Figure 3 plots K(S) with respect to ?fid» 0 tb e ,- K value for the best fit 
Euclidean model, 1.48x10“^ (counts s”* cm -2 ) 1 . Other Euclidean models 

would appear on Figure 3 as horizontal lines, K(S) constant. The 

trumpet-shaped region on the right shows the behavior of power law 
models acceptable at the 90% level. The hatched line at the far left is 
where the power law models must be terminated to avoid exceeding the 
total sky flux (equation [3]). With our data, an acceptable power law 

model that is truncated between the hatched line and the left hand edge 

of the. trumpet shaped region is statistically indistinguishable from the 
model as continued to the hatched line. The right hand edge shows the 
limits of the HEAP 1 resolved source counts. 


The greatest constraint placed by the fluctuations is on sources 
roughly an order of magnitude in flux bolow resolved sources. However 
care is required in interpreting the limits on N(S). The formal 
validity of the confidence region rests on the assumption that the 
actual N(S) Is well modelled by a single power law continuing without 
change of index past the lower limits of the region. More complicated 
models that do not lie wholly within the indicated region may be 
acceptable, J 2 .&. the dashed line of Figure 3. This line is a schematic 
representation of the N(S) behavior of sources observed directly at 
higher fluxes and extrapolated without evolution to the lower values of 
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S using the appropriate luminosity functions* 

To assess the degree of source evolution at low fluxes, direct 
measuments of N(S) would be the easiest to interpret. The deep surveys 
performed with the Einstein observatory [.e.^. Glacconi jet nl. 1979] 
provide such information. However Information from the deep survey, 
indicated by the two bars labeled DS in Figure 3, also have problems of 
interpretation. One difficulty is in transforming source fluxes from 
the Einstein 1-3 keV band to the higher energy band measured in the A-2 
data without good spectral information. The upper-right bar of the pair 
uses the published presumption that the sources have a 1.4 index power 
law photon spectrum, typical of the unresolved sky flux in the 2-10 keV 
band. The lower-left bar instead assumes that the deep survey sources 
have the 1.7 index spectra characteristic of active galactic nuclei. 
Both bars assume that N(S) is Euclidean at the deep survey flux limit. 
If y were nearer that of the ur evolved population's models at that flux, 
Y ~1.8, or if the sources were strongly evolving so that y were near 3, 
the bars would be adjusted to 50%-i30% of their indicated value. 
Conclusions drawn from the deep survey results must explictly consider 
the impact of these assumptions. Results from the Einstein medium 
survey [see £•£,• Maccacaro j^t jal. 1982] show source evolution less 
indirectly. 

6_. Excess Variance 

The fluctuations can be totally described by models of non- 
evolving sources, such as the dashed line in Figure 3. The 90% upper 
bound to any additional variance, added as a pure Gaussian,, is o 2 j S 
0.057 (counts exp” 1 ) 2 . If we assume that the unevolved populations 
account for 20% of the sky intensity then cj is SI. 7% of the '.remaining 
intensity. Any other source of variation is constrained by this limit. 

Any evolving populations 's distribution, N ev (S), is a source of 
fluctuations. It must satisfy the integral constraint 

a 2 s “ / dS S 2 N ev (S) . [5] 

c 2 s is a measure of the excess variance that is independent of the 
measurement solid angle, a 2 g £ 7x10**^ (counts s cm ) 2 . The impaCi of 
this limit depends on the particular form of N ev (S) but if the evolved 
sources are to make up the remainder of the sky flux we place a lower 
limit on the number of sources at 75 per square ^degree and an upper 
limit on their mean flux at 1.5x10 counts s cm (2.0x10 erg s 
cm” 2 ). 1^ we assume that N (S) is of the form 4 jt K qv S“ j we can set a 

limit K<1.7xl0“ 5 (counts s" t 'cm” 2 ) 2 , indicated by the right-hand hatched 
line in Figure 3. A wide .Latitude for the behavior of N ey (S) is 

allowed. 

These limits all assume that the sources that make up the 
background are not clustered, that is, their distribution among the 
measurements is Poisson. We can estimate that the total allowed 
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variation on a scale of 26 square degrees, roughly the size of a Schmidt 
survey plate, is 2.3*. If ^SOs contribute oil of the rematning^sky 
flux, then this limit means that the total (6 n/NWq $ 0.023, whore N is 
the moan number of QSOs on n oeale of 26 square negroes , and <$N is the 
total variation Including Poisson statistics and clumping. If we assume 
the number of QSOs is ^200 per square degree in order to estimate the 
Poisson noise portion, the additional variation due to clumping can be 
at most 1.9^. If QSOs are observed to have clustering with a larger 
value of on these scales, then the excess variance places an upper 
bound on their contribution to the total aky flux. 


Figure 4* Preliminary upper 
bounds on magnitude of large 
scale structure from limits of 
X-ray excess variance. 


Figure 4 shows estimates of the limits on the large scale 
structure of the universe placed by a previous upper bound on the excess 
variance. The three curves illustrate the dependence of the limits on 
models of the structure au well as on the origins of the X-ray sky. Two 
of the curves compare the case for an unevolved origin of the sky flux, 
z~0.5, and the density variation restricted to a single scale or defined 
by a power spectrum with index 1.8, typical of the inferred distribution 
at smaller scales [see je.jg,. Peebles 19801 « The third curve is for a 
single scale of variation where the sky flux is dominated by an evolved 
component with 2$zS3. 

JL* Conclusions 

All small scale fluctuations are consistent with what we expect 
from known, unevolved populations of sources resolved by HEAP 1. The 
large scale variation, dominated by i>\ galactic anisotropy, also admits a 
dipole signal which may be interpreted os a Gompt on-Get ting dipole 
signal consistent with the microwave results, although other 
interpretations are possible. The bound on any excess variance places 
limits on structures at moderate redshift and of large scale otherwise 
not easily accessible. 

We would like to thank our colleagues at GSFC, particularly 
C.M. Urry, for their contributions and critical advice on this 
presentation. ACF aknowledges the financial support of the Royal 
Society of London. 
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ABSTRACT 


HEAO I A2 and Einstein SSS spectral observations of Seyfert 
galaxies and BL Lac objects suggest; that in both cases, the X-ray 
emission is due to relativistic particles. The five BL Lac objects have 
very soft spectra and at higher energies (above 10 keV) may have hard 
tails. Combining our X-ray data with radio, infrared, optical, and 
ultraviolet observations, we can fit the BL Lac spectra with the 
familiar synchrotron self-Compton model if we allow for relativistic 
beaming (Urry and Kushotzky 1982, Urry et_ al . 1982). We show that 
Doppler beaming of an underlying (Seyfert-like) source population 
flattens the observed luminosity function, and we emphasize that the 
relative numbers of BL Lacs and quasars in given spectral intervals are 
strong functions of selection effects, the degree of Doppler beaming, 
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and the form of the Intrinsic luminosity function. 

The twenty-eight X-ray spectra of Seyfcrt galaxies are remarkably 
homogeneous: all are well-fit by power laws with mean energy spectral 
index << - 0.65 + 0.13, where the latter number indicates the dispersion 
(Mushotzky ejt al. 1980, Mushotzky 1982). This power law, taken with 
the IH-optlcal-UV emission, suggests that Compton scattering is the 
dominant X-ray production mechanism. No changes In spectral form are 
seen on short (10-10^ s) or long (0.5-1. 5 yr) timescales in either the 
A2 (2-40 keV) or SSS (0.5-3. 5 IceV) data even when the Intensity of the 
source changes. With one major exception (NGC 6814), no variability 
with AX/I > 0.1 is seen for timescales of 5 seconds to 6 hours (Tennant 
and Mushotzky 1982). For 6 month timescales, at most one third of the 
galaxies are variable in Intensity, and these tend to be the lower 
luminosity objects. 
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